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Abstract. We consider the mean time to absorption by an absorbing target of a 

diffusive particle with the addition of a process whereby the particle is reset to its initial 

position with rate r. We consider several generalisations of the model of M. R. Evans 

and S. N. Majumdar (2011), Diffusion with stochastic resetting, Phys. Rev. Lett. 106, 

160601: (i) a space dependent resetting rate r(x) ii) resetting to a random position z 

C^ ' drawn from a resetting distribution Viz) iii) a spatial distribution for the absorbing 

target Pt{x). As an example of (i) we show that the introduction of a non- resetting 

, window around the initial position can reduce the mean time to absorption provided 

^ i that the intial position is sufficiently far from the target. We address the problem of 

O \ optimal resetting, that is, minimising the mean time to absorption for a given target 

I ^i ' distribution. For an exponentially decaying target distribution centred at the origin 

we show that a transition in the optimal resetting distribution occurs as the target 

\N ' distribution narrows. 

> 

in 

(N 
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'^ . 1. Introduction 

l> 

Search problems occur in a variety of contexts: from animal foraging [1] to the target 

search of proteins on DNA molecules [2H1] ; from internet search algorithms to the more 

mundane matter of locating one's mislaid possessions. Often search strategies involve 

/\ ' a mixture of local steps and long-range moves |5H9]. For human searchers at least, a 

c^ ■ natural tendency is to return to the starting point of the search after the length of time 

spent searching becomes excessive. 

In a recent paper [10] we modelled such a strategy as a diffusion process with an 
additional rate of resetting to the starting point xq with rate r. Considering the object 
of the search to be an absorbing target at the origin, the duration of the search becomes 
the time for the diffusing particle to reach the origin. Statistics such as the mean time 
to absorption of the process then give a measure of the efficiency of the search strategy, 
defined by the resetting rate r. Moreover, the model provides a system where the 
statistics of absorption times can be computed exactly. 

A related model, where searchers have some probabilistic lifetime after which 
another searcher will be sent out, has been studied by Gelenbe [H] and mean times 
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to absorption computed. Also, in the mathematical literature the mean first passage 
time for random walkers that have the option of restarting at the initial position has 
been considered [T2] . 

In [To] it was shown that the mean first passage time (MFPT) to the origin for 
a single diffusive searcher becomes finite in the presence of resetting (in contrast to a 
purely diffusive search where the MFPT diverges). Moreover the MFPT has a minimum 
value as a function of the resetting rate r to the fixed initial position xq. Thus, there is 
an optimal resetting rate r as a function of the distance to the target xq. 

In this work we address the question of resetting strategies which optimise the 
MFPT in a wider context. To this end, we make several generalisations of single- 
particle diffusion with resetting studied in [10]. First, we consider a space dependent 
resetting rate r{x). Second, we consider resetting to a random position z (rather than 
a fixed xq) drawn from a resetting distribution V{z). Finally, we consider a probability 
distribution for the absorbing target Pt{x). The general question we ask is: what are 
the optimal functions r{x), V{x) that minimise the MFPT for a given Pt{x)7 Although 
we do not propose a general solution, the examples we study turn up some surprising 
results and illustrate that answers to the problem may be non trivial. 

The paper is organised as follows. In section 2 we review the calculation of the 
mean first passage time for one-dimensional duffusion in the presence of resetting to the 
intial position with rate r. In section 3 we introduce spatial dependent resetting r{x) 
and work out the example of a non resetting window of width a around the intial point. 
In section 4 we consider the generalisation to a resetting distribution V{z) and to a 
distribution of the target site Prix). In section 5 we formulate the general problem of 
optimising the mean first passage time with respect to the resetting distribution V{z). 
We consider the example of an exponential target distribution and show that there is a 
transition in the optimal resetting distribution. We conclude in section 6. 

2. First passage time for single particle diffusion with resetting 

We begin by briefly reviewing the one-dimensional case of diffusion with resetting to 
the initial position xq (see Fig. [1]), introduced in Ref. [10]. The Master equation for 
p{x,t\xo), the probability distribution for the particle at time t having started from 
initial position xq, reads 

dp{x,t\xo) (Pp{x^\xo) / H ^ , ./ N n\ 
oT = D — rp{x,t\xo)+rd{x-Xo) (1) 

with initial condition p{x,0\xq) = 6{x — Xq). In Eq. ([1]) D is the diffusion constant of 
the particle and r is the resetting rate to the initial position xq. The second term on the 
right hand side (rhs) of Eq. [1] denotes the loss of probability from the position x due to 
reset to the initial position xq, while the third term denotes the gain of probability at 
Xq due to resetting from all other positions. 
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Figure 1. Schematic space-time trajectory of a one dimensional Brownian motion 
that starts at xq and resets stochastically to its initial position xq at rate r. 



The stationary state of ([T]) is the solution of 

d^p*{x\xo) 



D- 



dx"^ 



rp*{x\xo) = —r6{x — Xq) 



(2) 



which is determined by the elementary Green function technique, which we now recall. 
The solutions to the homogeneous counterpart of ([2]) are e^"""^ where 

«o = V^ . (3) 

The solution to (|2]) is constructed from linear combinations of these solutions which 
satisfy the following boundary conditions: p* — )■ as a; — ?■ ±oo, and p* is continuous at 
X = Xq. Imposing these conditions yields 

p*{x\xo) = Aexp{—ao\x — xo\) ■ (4) 

Note that (jlj) has a cusp at x = xq. The constant A is fixed by the discontinuity of 
the first derivative at x = Xq which is determined by integrating ([2]) over a small region 
about xo 



9j9*(x|xo) 



dx 



dp*{x\xo) 



dx 



-an 



(5) 



p (x|xo) = Y^xpl 



Carrying this out yields A = ao/2 so that 

-ao|a;-xo|) • (6) 

Alternatively, the constant A in (jl]) could be fixed by the normalisation of the probability 
distribution @. 

Note that (|6]) is a non-equilibrium stationary state by which it is meant that there 
is circulation of probability even in the one-dimensional geometry. At all points x there 
is always a diffusive flux of probability in the direction away from xq given by —Ddp/dx, 
and a nonlocal resetting flux in the opposite direction from all points x 7^ xq to xq. 
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2.1. Mean first passage time 

We now consider the mean first passage time for tlie diffusing particle to reach the origin. 
One can think of an absorbing target at the origin which instantaneously absorbs the 
particle (see e.g. [T3]). 

A standard approach to first-passage problems is to use the backward Master 
equation where one treats the initial position as a variable (for a review see Ref. [E]). Let 
Q{x,t) denote the survival probability of the particle up to time t (i.e. the probability 
that the particle has not visited the origin up to time t) starting from the initial position 
X. The boundary and initial conditions are Q(0,t) = 0, Q(x,0) = 1 (see e.g. [15] for 
more general reaction boundary conditions). 

The backward Master equation (where the variable x is now the initial position) 
reads for the survival probability Q{x,t) 

dQ{x,t) d'^Q{x,t) 

— Qj. — = ^ — Q^ rQ{x,t)+rQ{xo,t) . (7) 

Note that Q{x,t) depends implicitly on the resetting position xq due to the third term 
on the right hand side of (jTj). The second and third terms on the rhs correspond to the 
resetting of the initial position from x to Xq, which implies a loss of probability from 
Q{x, t) and a gain of probability to Q{xo, t). 

Equation (jTj) may be derived as follows. We consider the survival probability 
(5(x, t + At) up to time t + At, where At is a small interval of time. We divide the time 
interval [0, t + At] into two intervals: [0, At] and [t, t + At]. In the first interval [0, At], 
there are two possibilities: (i) with probability rAt, the particle may be reset to Xq and 
then for the subsequent interval [At, t + At] this xq will be the new starting position 
and (ii) with probability (1 — rAt), no resetting takes place, but instead the particle 
diffuses to a new position (x + S,) in time At where ^ is a random variable distributed 
according to a gaussian distribution P{C) = (47rDAt)~^/^exp(— ^^/4DAt). This new 
position (x + ^) then becomes the starting position for the subsequent second interval 
[At, t + At]. One then sums over all possible values of ^ drawn from P(^). Note that we 
are implicitly using the Markov property of the process whereby for the second interval 
[At, t + At], only the end position of the first interval [0, At] matters. Taking into acount 
these two possibilities, one then gets 

Q(x, t + At) = r Atg(xo, t) + (1 - r At) j diP{i)Q{x + e, t) (8) 

which can be rewritten as 

Q(x,t + At)-Q(x,t) _ /d^p(^)^g^^^^^^)_g^^^^))^^g^^^^^)_^g(^^^)^Q(^^) ^^g) 



At J At 

Taking the limit At -^ then yields ([7]). 

The mean first passage time T to the origin beginning from position x is obtained 
by noting that — g^' ' dt is the probability of absorption by the target in time t — )■ t+dt. 
Therefore, on integrating by parts, we have 



/' 



CO ri( )iT f\ roo 



t^%^dt= / g(x,t)dt (10) 

ot Jo 
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D^^-rT{x)+rT{xo) (11) 



(assuming that tQ{x,t) — )• as t — )■ oo). Integrating ([7]) with respect to time yields 

dx"^ 
with boundary conditions T(0) = and T{x) finite as x — )■ oo. 

To solve for the mean first passage time beginning at the resetting position x = Xq 
we first consider the initial position to be at x > 0, different from the resetting position 
Xo, then solve ( ITT]) with arbitrary x and xq. Once we have this solution we set x = xq 
to determine T(xo) self-consistently. 

The general solution to (ITTl) is 



T{x) = Ae"o" + Se-"«" + ^ ^ ""^^^"^ (12) 



where Uq = Jr/D. The boundary condition that T{x) is finite as x — >■ oo implies A = 
and the boundary condition T(0) =0 fixes B. Thus 

T(x) = ^^""^^^"^ [l - e-"»i . (13) 

r L -I 

Solving for T(xo) self-consistently yields 



T(xo) = - [exp(aoa;o) - 1] = - 



exp 



\A7^^o) - 1 



(14) 



Note from (fT4ll that, for fixed Xq, T is finite for < r < oo. As a function of r for 
fixed Xo, T diverges when r — )■ as 

This is expected since as r — )■ 0, one should recover the pure diffusive behaviour (no 
resetting) for which the T is divergent-due to the large excursions that take the diffusing 
particle away from the target at the origin. Also T diverges rapidly as r — )■ oo, the 
explanation being that as the reset rate increases the diffusing particle has less time 
between resets to reach the origin. In other words, the high resetting rate to xq cuts off 
the trajectories that bring the diffusing particle towards the target. 

We now consider T as a function of r for a given value of xq and define the reduced 
variable 

/ ^ \ 1/2 
Z = aoXo = i—j Xq. (16) 

Since T diverges as r — )■ and r — t- oo it is clear that there must be a minimum of T 

irri 

with respect to r (see Fig. [2]). The condition for the minimum, — — = 0, reduces to the 

dr 
transcendental equation 

^ = l-e-^ (17) 

which has a unique non-zero solution z* = 1.59362.... In terms of the restting rate, 
this means an optimal resetting rate r* = {z*)'^D/xq = (2.53962 . . .)D/xq, for which the 
mean first passage time T(xo) is minimum. The dimensionless variable z (TT6l) is a ratio 
of two lengths: xq, the distance from the resetting point to the target, and (D/r)^/^, 
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Figure 2. The mean first passage time T — ^ cxp I \/r/D xq J — 1 plotted as a 
function of r for fixed xq = 1 and D ~ 1. clearly T diverges as r — >■ and as r — >■ oo 
with a single minimum at r* = 2.53962 . . .. 



which is the typical distance diffused between resetting events. Thus, for fixed D and 
Xq the mean first passage time of the particle can be minimised by choosing r so that 
this ratio takes the value z*. 



3. Space-dependent resetting rate 

In this section we generalise the model of section 2 to the case of space-dependent 
resetting rate r{x). 

The master equation for the probability distribution p{x,t\xo) is generalised from 
©to 

dp{x,t\xo) ^d^p{x,t\xo) , . , I N ,/".// /w ' +1 \xf \fM 

— —r[x)p[x,t\xo)+ ax r[x )p{x ,t\xo) d{x—Xo)[li 



dt dx^ 

The third term on the right hand side now represents the flux of probability injected at 
Xo through resetting from all points x ^ Xq. 

The stationary distribution p* (a; |a;o) satisfies 

d'^p*{x\xo) 



D- 



dx"^ 



r{x)p*{x\xQ] 



dx' r{x')p* {x'\xq) 5{x — Xq) 



(19) 



In general the stationary state is difficult to determine unless r{x) has some simple 
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D^^ - r{x)T{x) + r(x)r(xo) (20) 



form. The equation for the mean first passage time becomes 

d^T{x) 
dx"^ 
Again, this is difficult to solve generally for arbitrary r(x). 

In the following we consider a solvable example where r{x) is zero in a window 
around Xq and is constant outside this window. 

3.1. Example of a non-resetting window 

We consider the case of a non-resetting window of width a about |a;o|, within which the 
resetting process does not occur: 

r(x) = for |x — Xo|<a (21) 

= r for |x — xo| > a . (22) 

This choice is a rather natural one in the sense that a typical searcher usually doesn't 
reset when it is close to its starting point, but rather the resetting event occurs when it 
diffuses a certain threshold distance a away from its initial position. 
The Master equation reads 

^D ^-^ V rh[t)d[x — xq) \x — xq\ < a (23) 



dt dx 

dx'^ 



«9^p(a;, tlxo) / i n i i 

D — rp(x,t|xo) |x-Xo|>a (24) 



where 

h{t) = / dxp(x, t|xo)6'(|x — Xol — a) , (25) 



with initial condition p(x, 0) = 5{x — Xq). Thus h{t) is the probability that the particle 
is outside the non- resetting window, i.e., in the resetting zone at time t; the particle is 
reset to the origin with a total rate h{t)r. 

First, we consider the stationary state. One can solve for the stationary probability 
using the Green function technique of section 1. For |a; — xo| > a (outside the window), 
p*{x\xo) satisfies D — ''^ q^2 = rp*{x\xQ) and should tend to zero as |x| — t- oo. For 
< |x — xo| < a (inside the window), p*{x\xq) satisfies D — ^q^2 = for all x ^ Xq. 
The solution should be continuous at x = Xq, but its derivative must undergo a jump 
at X = Xq and the jump discontinuity can be computed by integrating Eq. |23] across 

X = Xq. 

Thus, noting that the solution should be symmetric about x = xq, one has 

p*(x|xo) = Aexp — ao(|x — xo| — a) for |x — xo|>a (26) 

= y4 — i? (|x — Xol — a)) for |x — Xo|<a (27) 



where ao = Jf/D and the constants A and B are determined by the discontinuity in 
the derivative of p*(x|xo) at x = xq and the continuity of the derivative at |x — xo| = a. 
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Figure 3. The stationary solution p*(x|a:o) in Eq. EM?7l plotted as a function of x, for 
the choice a;o = l,a = l,r = l and D = 1. The nonresetting window is over x € [0,2] 
with the initial position at xq = 1. The solution is symmetric around xq = 1 with a 
cusp at X = xo = 1. 



The result is 
A = 



an 



B = aoA 



(2^ 



2 + 2aao + a^al 

The solution has a cusp at x = xq and a discontinuity in the second derivative at 
|x — xo| = a (see Fig. [3]). 

We now consider an absorbing trap at the origin. The backward equation for T(x), 
the mean time to absorption beginning from x, reads 

d^T(x) 



D- 



-1 = D 



dx'^ 
d^Tjx) 

dx'^ 



rT{x)+rT{xo) for |x — Xo|>a 



for Ix — Xo| < a 



The general solution to fl29|30p is 



T{x) = Ax + B 



X 



2D 



and the solution that does not diverge as x — )■ oo is 

1 + rT{xo) 



T{x) 



_|_ /~i -ao{x—xo-a) 



^ 1 + rT(xo) ^ ^^_„„, ^ ^^„„, 



for X > xo + a 
for X < xo — a 



(29) 
(30) 

(31) 

(32) 
(33) 
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The constants A, B, C, E, F are determined by the continuity of T{x) and T'(x) at 
|x — XqI = a and the boundary condition T(0) = 0. The result for T{xq) is 

1 



T(xo) 



r(l + aao) 



coshao(3;o — a) \1 + Zaa^ -\ h a a^ 



+ sinhao(xo — a) 1 + 2aao 



Sa^tto' 



(34) 



2 

We now consider the reduced parameters z = OoXq and y = a^a, and T as a function 
y for z fixed. The allowed values of y are < y < z. At y = z, one can show that 
> 0. Therefore the minimum of T with respect to y is either at y = or at a 



dT 

dj/ 



i/=z 



non trivial minimum < y < z. 

The condition for a minimum 4^ = reduces to 

ay 



tanh(2; — y) (35) 



4: + 5y + 2?/2 

Therefore the condition for there to be a nontrivial minimum for y > is given by 
tanh^ > 1/2 or equivalently z > (log3)/2 = 0.5493 . . .. 

In summary, the analysis of the condition for T{y) to be a minimum reveals that: 
if 2; < (log3)/2 then y = is the minimum of T{y)] if z > (log3)/2 then T{y) has a 
nontrivial minimum at < y < z. Therefore, when z < (log3)/2 the introduction of 
a window around the initial site where resetting does not take place does not reduce 
the mean time to absorption. A strategy of introducing a non-resetting window is 
an effective one only when the initial point is sufficiently far from the search target. 
Otherwise it is advantageous to always reset. 

3.2. Optimal resetting function 

Having seen in the previous example that non-trivial behaviour emerges for a simple 
spatial-dependent resetting rate r{x), one can ask for the optimal function r(x). The 
optimisation problem would be to minimise T under certain constraints pertaining to 
the information available to the searcher. Clearly if there are no constraints, that is 
one can use full information about the target position, the optimal strategy is to reset 
immediately whenever x > xq and not reset when x < xq. This corresponds to the 
choice 

r(x) = for X < Xq 
r{x) = 00 for X > Xq . 

In this case problem ( l20l) reduces to the mean first passage time of a diffusive particle 
with reflecting barrier at xq the solution of which is 

T*{xo) = ^ . (36) 

Thus, f l5B]l gives the lowest possible mean flrst passage time for a diffusive process. One 
can then ask about how close simple strategies, such as a spatially constant resetting 
rate r or non-resetting window, come to approaching this bound. 
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For example, the case of spatially constant resetting rate r considered in section 2 
yields a minimum MFPT using ([1] 



T = |1^ = ^- = 3-0883... T*(xo) (37) 

As noted in section 3.1 the value 3.0883 may be improved upon by considering a non- 
resetting window around Xq. 

However, (l36ll uses the crucial information of whether the target (at a; = 0) is to 
the right or left of the resetting site Xq. More realistically, the searcher would not have 
this information. The relevant optimisation problem is to find the optimal resetting rate 
r(|x — xo|) (constrained to be a function of the distance |x — xo| from the resetting site) 
that minimises T{xq). This remains an open problem. 

4. Resetting distribution and target distribution 

In this section we consider the generalisation to a system with resetting to points 
distributed according to V{z). We shall also consider a distribution of the target site 
Pt{x). 

4-.1. Stationary state 

We begin by considering again the one-dimensional case of diffusion but this time with 
resetting to a random position: at rate r the particle is reset to a random position 
z -^ z + dz drawn with probability V{z)dz. We refer to V{z) as the reset distribution. 
For simplicity we take the initial position xq to be distributed according to the same 
distribution as the reset position p{xo, 0) = V{xo). 

The Master equation for the probability density p{x, t) now reads 

dp{x,t) d'^p{x,t) 

— Q-^ — = D^-^ rp{x,t) + rV{x) . (38) 

The stationary solution to fl38|) is simply found using ([6]) as the Green function: 

p*(x) = f dzV{z)p*{x\z) (39) 

which, using p*{x\xo) given by ([6]), yields 

p*(x) = Y [ dzViz)exp{-ao\x - z\)) . (40) 

4-2. Mean first passage time 

The mean first passage time, T{xo,xt), to a target point xt, starting from xq with 
resetting distribution V{z), satisfies 

d^T{xo,XT) 



D- 



dxl 



rT{xo, Xt) + r dzV{z) T{z, xt) (41) 

with boundary condition T{xt, xt) = 0. To solve this equation we let 

F{xt) = / d2 V{z)T{z, Xt) (42) 
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then write down the general solution to (141 p and solve for F{xt) self-consistently. 
The general solution of (HTj) which is finite as Xq — > c>o is 



T{xo,xt) = Ae-"«l^°-^^l + - + F{ 



Xt) 



dzV{z)e 



The boundary condition T{xt, xt) = implies A 
expression for A in (143!) and integrating we find 

F{xt) =(- + F{xt 
\r 

which yields 

FM = - f ^^ 

r \2p*[XT) 

Inserting this into ( 143|) we obtain 



-aa\z—XT\ 



T{xo,xt) 



"0 



[1 — exp(— ttol^o ~ xt\ 



(43) 
Then substituting this 

(44) 
(45) 
(46) 



2rp*{xT) 

As noted above it is convenient to choose the same distribution for xq as the resetting 
distribution. Averaging over xq then gives using (!39l) 



T{xt) 



ttO 



(47) 



2p*{xt) 

Equation ( 147|) gives the expression for the mean first passage time to a target positioned 
at Xt- Let us check the case of a single position xq to which the particle is reset 
V{z) = 6{z — Xq). In this case (H71) becomes 



T{xt) = - 
r 



«o 



1 



1 



(4J 



2 p*{xt) 

which recovers (IT^ when xt is set to 0. 

Finally, we average over possible target positions drawn from a target distribution: 
Pt{xt) 



r 



'^fdxr^4^-l 
2 J p*[xt) 



(49) 



Equation (H9|) gives the main result of this section — the MFPT for a resetting 
distribution V{xq) and averaged over target distribution Pt{xt)- 



5. Extremisation of mean first passage time 

Let us now consider the problem of extremising T given by (l49l) . for a given target 
distribution Pt{x), with respect to the resetting distribution V{z). Throughout 
this section we will assume a symmetric target distribution: Pt{x) = Pt{—x) and 
P^{x) = -PH-x). 



The problem is to minimize the functional appearing in ( H9l) : / da: L^? where 



P*{x) = -^ dzV{z)e~ 



■ao\z—x\ 



(50) 
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subject to the constraint J dzV{z) = 1. The functional derivative to be satisfied is 

(51) 



(52) 



6 
SV{y) 


Ux^^[''}+\ [dxV{x) 
J p*{x) J 


= 


where A is a Lagrange multipher. Condition (51 1) yields 


J [p*{x)] ao 


For (I52l) to hold for all y requires that 


Pt{x) 


- = A, 


p*{x)f 


or fixing A through the normalisation oi p*{x) 


p\x) = 


P^^) 




JdzP^/\z) ■ 



(53) 



(54) 

Equation (IMI) implies that to minimise T the stationary probability distribution 
should be given by the square root of the target distribution. This result has been derived 
in [16] for the case of searching for the target by sampling a probability distribution V{x). 
This corresponds to the limit r — ?> oo of our model. For r < oo we have the additional 
constraint that the optimal p* should be realisable from a resetting distribution V{z) 
through formula ( !50l) . 

Equation ( l50l) may be solved for V{z) for a desired p*{x) by taking the Fourier 
transform and using the convolution theorem to give 

V{k)=(l + ^]p*ik) (55) 

where V{k) is the Fourier transform of V{x) and p*{k) is the Fourier transform oi p*{x). 
We may invert the Fourier transformation to find 

V{x) =p [x)- — . 56 

Oq dx^ 

However this solution may become negative in which case the solution to the 
optimisation problem is unphysical. 

5.1. Example of an exponential target distribution 

As a simple example, we consider an exponentially decaying target distribution peaked 
at X = 0: 

PAx) = ^e-^l^l . (57) 

We first note that for a delta function resetting distribution V{z) = 6{z — xq) the mean 
first passage time (1151) diverges when oq > /?. Therefore, for small /3 (a broad target 
distribution) one expects an optimal resetting distribution (for fixed ao) that differs 
from a delta function. 
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For /3 < 2ao, the optimal stationary distribution is from fl54p 

p*^^) = ^e-^l^l/2 (58) 

This expression yields from fl56l) a resetting distribution that is always positive, thus 
the optimal resetting distribution 



V{z) = ^e 



-/3|2|/ 



4a^ 



i|*(.) (59) 



'0 

For /3 > 2ao, f l59|) always gives negative probabilities due to the first term. Therefore 
we anticipate that V{x) = 6{x) is at least a locally optimal solution. In fact one can prove 
this is the case by showing that any distribution of the form V{x) = (1 — e)(5(x) + ef{x), 
where f{x)>0 and Jdxf{x) = 1 leads to an increase in (H9i) at first order in e when 
P > 2aQ. (As the proof is straightforward but somewhat tedious we did not include it 
here.) Thus a transition in the form of the optimal resetting distribution, from a single 
delta function to ( l59l) . occurs at /3 = 2ao. 



5.2. Inversion ofp*{x) 

As noted above, the constraint V{x) > means that the optimal p*{x) given by (l54l) 
may not be realisable from a physical resetting distribution V{z). We are therefore led 
to the general question of when a desired stationary distribution (e.g. ( !54l) ) which we 
denote g{x) may be generated from (150|) i.e. when can we invert 



g{x) = yJ d2P(2)e-°°l^-"l (60) 

to obtain a physical V{z)l 

Let us first discuss a sufficient condition for the resetting distribution implied by 
fl60|) to be physical. 

Equation ( l55l) relates the characteristic functions of the two distributions V{x) and 
g{x) (given there by p*(x)). In terms of the characteristic function, Polya's theorem [17] 
states that if a function 0(/c) satisfies: 0(0) = 1; </)(fc) is even; </)(A;) is convex for /c > 0, 
and 0(oo) = 0; then (/)(/c) is the characteristic function of an absolutely continuous 
symmetric distribution. Polya's theorem therefore gives a sufficient condition for V{x) 
implied by p*{x) to be physical. 

The condition for convexity becomes in one dimension 

' "^ >0 for all k>Q. (61) 



dfc2 



«n 



9{k) 



If the function V{k) does not satisfy the conditions of Polya's theorem, the solution of 
fl60|l is invalid as a probability distribution i.e. the desired g{x) cannot be realised from 
any resetting probability distribution V{z). 

In the case where (160|) may not be inverted to give a physical V{x), it may be 
possible to generate the desired form for g{x) on a finite region by choosing a compact 



Diffusion with Optimal Resetting 



14 



support for V{z). Let us assume g{x) to be a symmetric function of x. Then if we 
choose 

1 d'giz) 



9{z} 



Oq dz"^ 



for \z\ < yo 



V{z) = X 

= for \z\ > yo , 
where A is a normahsing constant, we find 
p*{x) = \g{x) for |x| < yo 
p*{x) = A^(i/o)e"''(^°"l^l^ for |x| > yo 
provided that yo is chosen so that 

9{yo) + —g'iyo) = o 
9{-yo) g'i-yo) = o . 

«0 



(62) 
(63) 

(64) 
(65) 

(66) 
(67) 



see 



Appendix A] ). The second condition follows from the first by the assumed symmetry 

(68) 



of g{x). As an example, we consider the gaussian distribution 

1/2 

g[x 



I; ^ 



./3x2 



The inversion of fl^ using fISBl) yields 

1 d'gix) 



V{x) = g{x) - 
which becomes negative for 



tto dx'^ 



'^S'\-^- 



n 



2/3 4/32^2 



\x\ > 



«0 

2/3 



1 + 



2/3^ 



On 



1/2 



However, choosing a compact support for P(^;) according to fl67|) . yields 

^°=2^ 
and we find that the resulting distribution (170|) is positive for all x. 



(69) 



(70) 



(71) 



6. Conclusion 

In this paper we have considered some generalisations of diffusion with stochastic 
resetting to the case of spatial-dependent resetting rate and a resetting distribution. 
We have considered the mean first passage time to a target which may be situated at a 
fixed point (the origin) or distributed according to a distribution and derived the result 
(H9|) . The minimisation of this quantity may then be formulated as an optimisation 
problem of which we have studied some examples. 

In particular we have seen some perhaps unexpected results. First, the introduction 
of a non-resetting window around a fixed resetting position reduces the MFPT when the 
target is sufficiently far away. This suggests that the optimal resetting distribution, in 
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the case where we consider a resetting rate that is symmetric about the restting point, 
r(|x — XqI) may be non-trivial. We have also seen that in the case of an exponentially 
distributed target ( 1571) the optimal resetting distribution undergoes a transition from 
(!59|) to a pure delta function at the origin. 

Generally, the computation of an optimal resetting distribution is an open problem 
since the resetting distribution that minimises T may be become negative over some 
domain and therefore nonphysical. In the case where (1561) becomes unphysical, although 
we do not have a solution to the extremisation problem of minimising T subject to the 
additional constraint V{x) > we may propose likely candidates for extremal solutions. 
One possibility for the optimal physical solution is one that has compact support i.e. 
since the constraint for the distribution to be physical is that V{x) > 0, one might 
expect that the optimal solution lies on the boundary where V{x) =0 for some regions 
of X. However, we have no proof that this is the case. 

Further considerations for optimising mean first passage times in a more realistic 
search process would be to add a cost to resetting since in the present model the 
diffusive particle instantaneously resets to its selected resetting position. This could 
be implemented by attributing some time penalty to each resetting event, as is the case 
in the framework intermittent searching. 

Appendix A. Proof that ([63]) yields ( 165|) 



We wish to show that expression (!63|) for V{x 

1 d'g{z 



r{z) = X 



9[z) 



otherwise , 



for \z\ < yo 



(A.l) 
(A.2) 



yields ( l64|) - (!65!) for the stationary distribution given by (140|) . provided that (l66l) - (!671) 
holds. 

We begin by inserting flA.ip - flA.2|) into ( HOjl in the case |a;| < yo: 



p*{x) 



9{z) 




dz2 



On 



-ao{x-z) 



9[z)- 



1 d'g{z) 



-ao{z-x) 



We use the following integration by parts, valid for all ao 7^ 



9{z) 



l_fg{z) 
tto dz^ 



'dz 



9{h) 



1 dg{z) 



ao dz 



z=b. 



^aoh 



"0 



9W 



1 dg{z) 



ao dz 



Inserting this into (1A.3P and cancelling terms yields 

-ao{x+yo) 



p [X) 



«oA J 2 g{x) e- 



2 I ao 

gao(a;-J/o) 

+ 

"0 



CtQ 

-9{yo) - 



9{-yo) 



1 dg{z) 



1 dg{z) 



ao dz 



-yo. 



tto dz 



z=yo. 



(A3) 



3 00a 



"0 



-(A.4) 



(A.5) 
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Then conditions (1M]) -( 1^ ensure that p*{x) = Xg{x) for |x| < yo- 
In the case x > yo we find 



p*{x) 



QfoA 



yo 



9{^) 



^Q^oyo 



-002; 



1 ^'gjz) 



am) 



-aQ(x-z)^^ 



ao dz 



z=yo. 



2 I ao 

= Ae-"«(^-^«)^(i/o) 

where conditions ( 166|) - (|67|) have been used. 

Similarly in the case x < —yo we obtain p*{x) = \g{—yo)e' 



-aoyo 



ao 



9[-yo) 



1 dgjz) 
ao dz 



ao(yo+x) 



-yo_ 



(A.6) 
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